Material-specific gap function in the high-temperature superconductors 
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We present theoretical arguments and experimental support for the idea that high-T c superconduc- 
tivity can occur with s-wave, d-wave, or mixed-wave pairing in the context of a magnetic mechanism. 
The size and shape of the gap is different for different materials. The theoretical arguments are based 
on the t- J model as derived from the Hubbard model so that it necessarily includes three-site terms. 
We argue that this should be the basic minimal model for high-T c systems. We analyze this model 
starting with the dilute limit which can be solved exactly, passing then to the Cooper problem which 
is numerically tractable, then ending with a mean field approach. It is found that the relative stabil- 
ity of s-wave and d-wave depends on the size and the shape of the Fermi surface. We identify three 
striking trends. First, materials with large next-nearest- neighbor hopping (such as YBa2Cu307_ x ) 
are nearly pure d-wave, whereas nearest-neighbor materials (such as L^-^Sr^CuO^ tend to be 
more s-wave-like. Second, low hole doping materials tend to be pure d-wave, but high hole dop- 
ing leads to s-wave. Finally, the optimum hole doping level increases as the next-nearest-neighbor 
hopping increases. We examine the experimental evidence and find support for this idea that gap 
function in the high-temperature superconductors is material-specific. 



PACS numbers: 74.20.De, 74.20.Mn, 74.72.-h 
I. INTRODUCTION 

In the years following the discovery of the high-T c ma- 
terials, theoretical debate often centered on the anoma- 
lous normal state properties, these being seen as the key 
to understanding the underlying physics. The supercon- 
ducting state, by contrast, was generally thought to be 
conventional except for its coupling strength. More re- 
cently, however, it has seemed more reasonable to regard 
the superconducting state, especially the gap symmetry, 
as holding the key to theoretical understanding. Recent 
experimental progress gives hope that the gap symmetry 
can be unraveled. Once this is done, one expects that 
very strong constraints can be placed on the microscopic 
model. This will certainly be true if, as we shall contend 
in this paper, different high-T c systems have different gap 
functions. 

The early consensus that the high-T c superconductors 
are s-wave has been replaced by the view that they are 
most likely d-wave. Support for this latter position has 
come primarily from three experimental sources: pene- 
tration depth measurements jl]|2|, photoemission mea- 
surements |^^], and Josephson interference measure- 
ments [§-0. These results very much strengthened the 
idea that the mechanism is magnetic. Even before the 
discovery of high-T c superconductivity it was clear that 
magnetic interactions based on antiferromagnetic corre- 
lations would quite naturally give rise to a tendency to- 
ward higher-wave pairing [^|,^), and studies of the t-J 
model [ fn)|]Ti"| and the Hubbard model [Q around the 
time of the discovery confirmed this for the CuC>2 square 
lattice . A more complete theory, though one which 
requires some phenomenological input, has been con- 
structed on the hypothesis that antiferromagnetic spin 



fluctuations act very much like the phonons in low-T c ma- 
terials and cause the pairing |i~3|]i~4| ]. This theory, which 
we shall call the spin-fluctuation model, requires the pres- 
ence of very strong spin correlations which lead to a high 
critical temperature and to d-wave pairing. 

If it is accepted that a magnetic mechanism is respon- 
sible for high-T c , there still remain many outstanding is- 
sues. The spin- fluctuation model generally relies on a 
pairing interaction which is very well localized in momen- 
tum space. The interaction is proportional to the mag- 
netic susceptibility, which is taken to be very strongly 
peaked near the k x = ±n/a and k y = ±7r/a points. 
The spin-fluctuation model leads unambiguously to d- 
wave superconductivity and not to s-wave. This also ap- 
pears to be the case in the spin bag model, although the 
range of the interaction is shorter in this case, of order 
hvp /Asdw, where vp is the Fermi velocity and Asdw 
is the charge density wave gap of the ordered phase [^5) . 
On the other hand, theories based on the t-J model rely 
on a spin-spin interaction which is rather local in real 
space. We shall call this the spin-interaction model. If 
the gap equation is solved in this model, substantial re- 
gions of both s-wave and d-wave pairing are found [fl0[ . 
In earlier variational Monte Carlo (VMC) calculations 
on this model, the dominant instability of the normal 
state was toward d-wave pairing, but in some parame- 
ter regimes there was also the possibility of admixture 
of s-wave symmetry into the ground state Jr|. Thus, 
although the spin-interaction model clearly falls into the 
category of magnetic mechanisms, it is distinct from the 
spin-fluctuation model in that it predicts that the gap 
symmetry is not necessarily always pure d-wave. Rather, 
s-wave and d-wave are competing instabilities. This lat- 
ter direction of inquiry has been briefly summarized by 
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Muller, who notes that the possible coexistence of s- and 
d-wa ve c an describe recent experimental results fl7j| . Our 
Sec. VII is an expansion of this theme. 



The comparison of the spin-fluctuation model and the 
spin-interaction model leads to the conclusion that all 
magnetic mechanisms are not alike. Even if it is conceded 
that the basic mechanism for high-T c superconductivity 
is magnetic, there is still work to be done. This paper is 
written to amplify and sharpen this point. 

We present a brief discussion of recent theoretical work 
which is directly relevant to the issue material-specific 
gap functions (Sec. |l]) arising from the magnetism mech- 
anism. The extended t-J model is developed from the 
overlapping copper and oxygen orbitals in the Cu02 
planes (Sec. III). We point out physical reasons to sup- 
port the idea of s- and d-wave as competing instabilities 
based upon the solution of the t-J model in the dilute 
limit (Sec. IV). It is noted that this solution in fact cor- 
responds to a physical limit of the gap equation in which 
the interactions do not possess a frequency cutoff. We 
then extend our arguments by solving the Cooper prob- 
lem with the t-J model (Sec. |y|). We follow with a pre- 
sentation of a mean field approach to the t-J model which 
does indeed yield mixed s- and d-wave pairing under cer- 
tain circumstances (Sec. VI). Finally, we survey existing 
experiments and find evidence of t he trends identified in 
the analysis of our model (Sec. VII). The approach of the 
experimental survey is to examine different materials one 
by one in order to see if they all have the same gap sym- 
metry. Our conclusion is that the details of the gap func- 
tion are material-specific, and that this lends support for 
using the t-J model, extended in an appropriate fashion, 
as a basic, yet flexible description of high-temperature 
superconductivity in the cuprate-layered materials. 



II. THEORETICAL BACKGROUND 

Admirably complete and comprehensive surveys of the- 
oretical work on the magnetic mechanism and the ques- 
tion of d-wave pairing has recently been carried out by 
Scalapino |18| (see particularly Appendix A, and refer- 
ences therein), Dagotto Jll|, and von der Linden [po[ , 
who stress the numerical work. They conclude that nu- 
merical evidence from quantum Monte Carlo, variational 
Monte Carlo, and exact diagonalization on small sys- 
tems favors d-wave pairing. Certain diagrammatic stud- 
ies confirm this. These references also point out that 
the evidence is far from conclusive. For example, quan- 
tum Monte Carlo searches for a finite superfluid density 
in the two-dimensional Hubbard model have not been 
successful. Thus, studies of d-wave superconductivity in 
magnetic models are highly suggestive, but we are not 
assured that superconductivity exists in these microscop- 
ically justified models. 

Our approach in this paper is to sidestep this thorny 
issue. Instead of investigating the microscopic mod- 



els with the most sophisticated tools available to ask 
whether superconductivity is present, we shall use sim- 
ple mean-field-type methods which are expected to give 
superconductivity. The idea is to apply these methods 
to rather more complicated models, intended to more 
closely mimic the actual systems. If we can identify 
trends as a function of the (fairly numerous) parameters 
in the models, these may then be compared to experi- 
ment in a detailed way. Accordingly, we do not review 
those studies which have been carried out on the nearest- 
neighbor Hubbard model or its strong-coupling equiva- 
lent, the nearest-neighbor t-J model, the aim of which 
has generally been to search for pure d-wave supercon- 
ductivity. Our interest is in the question of whether the 
gap structure may be material-specific. Specifically, we 
wish to address the issue of whether the shape and size 
of the gap function may depend on details of the band 
structure and the doping level. We therefore review here 
those papers which treat more complex models or more 
complicated gap functions. 

The finding of a more complex gap structure in the t-J 
model originates from the very beginning of this field of 
study. Ruckenstein et at, were among the first to apply 
this model for use in calculating some of physical parame- 
ters of the high-temperature superconductors jl(|. They 
analyzed the model with a mean field treatment (similar 
to Sec. VI), and found a parameter range where a mix- 
ture of s-wave and d-wave was the most stable ground 
state. 

Gros, et al, used a variational Monte Carlo method by 
which the ground state of the t-J model was studied with- 
out the need for further approximations such as mean 
field |IlJ]. This method has the advantage of explicitly 
handling the requirement of no double-occupancy. Fur- 
thermore, this method is a test for evaluating candidate 
wavefunctions and comparing their ground state energies 
with one another through extensive parameter surveys. 

Li and the present authors utilized this VMC method 
in studying the t-J model along with the three-site terms 
to suggest the possibility of sd-mixing in the cuprates 
ifjjjfl . In this paper, we presented VMC calculations 
which compared ground state energies of d-wave, ex- 
tended s-wave and various mixtures of the two. The 
results showed, as is repeated often throughout this pa- 
per, that d-wave was the ground state near half-filling, 
whereas an admixture of s and d won out at higher hole 
doping. This VMC work, done solely on tetragonal lat- 
tices, found little difference in energy among the various 
kinds of sd-mixing, namely s+d, s+id, or various phases 
in between. Rather, all such mixtures were deemed de- 
generate within the statistical error of the calculations, 
as can clearly be seen in Fig. 6(a) of that work. We 
proposed at the time that the introduction of anisotropy 
would lift this degeneracy, and indeed recent calculations 
indicate (with allowance for further necessary study) that 
in orthorhombic lattices, the VMC calculations favor an 
s+d- wave state [^T|. These are preliminary results of a 
nonsystematic study of the t-J model phase diagram, the 
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details of which are not clear at this time and demand 
further analysis before more definite claims can be made. 

The recent work of O 'Donovan and Carbotte yields s- 
and G?-wave mixing p2fl . Using a tight-binding model 
with orthorhombic distortion, they solved the zero- 
temperature BCS gap equation with a nearest-neighbor 
interaction, and found, as we do throughout this work, 
that only s x 2 +y 2 and d x 2_ y 2 components are found. In- 
deed, nearer to half-filling, they saw a predominant d- 
wave with small s-wave admixture, whereas at smaller 
density n, the s-wave component increased until it was 
the predominant phase with small d-wave admixture. 
The relative phases between the two was (a d+s-wave 
phase) shifting rather quickly as n decreased further to ir 
(s-d-wave phase) with the increase of the s-wave compo- 
nent. This cross-over of predominant d- to s-wave is in 
rough agreement with the trends noted in Jll| and also 
below. 

Beal-Monod and Maki also solved the zero-tempera- 
ture gap equation using a fermion-fermion interaction, 
where only states on the Fermi surface were considered 
in their analysis J23]. They found that even a small 
anisotropy enhanced the maximum gap as well as the 
transition temperature. In the tetragonal case, they 
found pure d-wave as the solution. 

Varelogiannis, in study of electron-phonon interactions 
in the 2D superconductors, has found that cither s- or 
d-wetve results from consideration of the Coulomb pseu- 
dopotential as well as doping pi] . That some electron- 
phonon interaction investigations are yielding similar 
qualitative results as spin-interaction models is interest- 
ing. 

The work most similar in spirit to the present is that 
of Dagotto and collaborators. They have also stressed 
that spin-interactions which are short-range in real space 
may lead to a rich superconducting phase diagram with 
a region of stability for s-wave |25| . These references 
also make the important point that details of the band 
structure, particularly peaks in the density of states, may 
have an important influence on superconductivity. The 
band structure dependence is also one of the important 
aspects of the results below; we find other, more local, 
aspects of the band structure to be important as well. 

Scalapino's survey indicates that most theoretical work 
on magnetic mechanisms has concentrated on the issue 
of "s-wave versus d-wave" : do the simplest models give 
rise to <i-wave superconductivity? The whole range of 
theoretical tools available has been applied to this prob- 
lem. If they do, and if the experiments demonstrate that 
the high-T c materials are <i-wave, then, taking everything 
together, we have evidence for a magnetic mechanism. 
Taking our cue from the papers which show more com- 
plex phase diagrams arising from magnetic interactions, 
we wish to add an additional element to this debate. Is 
it possible to identify trends in the theoretical results on 
the size and shape of the gap function which can be com- 
pared with experiment? Our strategy will be to look at 
the simplest magnetic model, the t-J model, but adding 



enough flexibility to the model that we can hope to make 
material-specific predictions. The model then becomes 
sufficiently complicated that we are limited to the sim- 
plest mean-field-type methods. These methods do have 
the advantages of simplicity and physical transparency. 




FIG. 1. Diagram of the CuO planes showing four unit 
cells and the overlay of the Cu 3d x 2 _ y 2 and O 2p x and 2p y 
orbitals. Note that a more realistic picture would show the 
actual hybrid orbitals of the high-T c CuO planes. 



III. THE MODEL 

In high-Tc systems, conduction takes place in the 
copper-oxygen planes. The unit cell consists of one cop- 
per and two oxygen atoms as seen in Fig. |l|. The chem- 
istry of these atoms suggest, and band structure calcula- 
tions confirm, that there is only one orbital on each site 
that plays a role in this conduction. These are the oxy- 
gen 2p x or 2p y orbitals (x or y depending on which one 
points to the neighboring copper atoms) and the copper 
id x 2_ y 2 orbital. The six states must accommodate five 
electrons in, for example, La2Cu04. We shall call this the 
half-filled band case, for reasons which will become clear 
later. In La2- x Sr x Cu04, they must accommodate 5 — x 
electrons, as Sr contributes 1 less electron per unit cell 
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to the CuO planes. The noninteracting model is easily 
solved. 

Let us take creation operators dj, p\ x , and p\ 2 - Then 
the Hamiltonian is 



-tpda^2 [( d lPrl +Pll d r) ~ (4Pr 
r 

+tpdaJ2 [( d UaxPvl +Pll d v+ax) 
r 

~( d r+ayPr2 + Pr2 d T+ay 



P\2 d v) 



(1) 



Here, Ed and e p are the atomic energies of the two orbitals 
and 



t 



pda 



d 3 r ip d (r) AH ip py (r - ay) 



(2) 



is the overlap integral of the 3d x 2_ y 2 orbital and the 2p y 
orbital in the unit cell. AH is the difference of the atomic 
and crystal Hamiltonians. All the other overlaps in our 
restricted basis set are related to this one by symme- 
try The relative signs in the Hamiltonian arc important, 
and are determined by the character of the atomic wave- 
functions. The above tight- binding Hamiltonian assumes 
that the orbitals are orthogonal to each other. 
The Bloch wavefunction is now defined by 



[«! pUk) + u 2 pt 2 (k) + u 3 4(k)] |0) , (3) 



with the itj(k) still to be determined. The matrix ele- 
ments of H are 



<v£iwivv> = 5^Y. M*) <w ^( k )' 



(4) 



with 
M k ) 



(5) 







W" 1 + e -* k v a ) t pda {l - e~^ a ) 



W" 1 + elkya ) 
t P da(l - e tk * a ) 
£d 



as may be verified by a simple calculation. The secular 
equation to determine the eigenvalues e(k) at a given 
wavevector is 



\h ij (k)-e(k)\=0, 
which takes the form 



(6) 



(s p - e(k)) 



(e p - e(k))(e d - e(k)) - Ul do sin 2 ' ° 



-4ip dff sin 2 



(7) 



There are three solutions to this equation, which we shall 
call £ n (k), e+(k), and £-(k). They are 



^n(k) — £p 3 

e„ + e d , 1 
2 



(8) 



(ed-e p ) 2 + m 2 pda sin 2 ( ^ 



16t pda . sin 



1/2 



(9) 



We see that the band corresponding to £„(k) is com- 
pletely flat. It has no amplitude on the Cu- atoms and 
can be thought of as a nonbonding orbital. £-(k) has 
dispersion which reaches downward from the atomic ener- 
gies and is a bonding orbital. e+(k) has dispersion which 
reaches upward from the atomic energies and is an anti- 
bonding orbital. It is generally found that e p w e d . This 
leads to the very strong hybridization between copper 
and oxygen which is characteristic of the high-T c materi- 
als and distinguishes them, chemically, from most other 
transition-metal oxides. The two energies lie about 3.2 
eV below the Fermi energy. (We shall not attempt to 
assign precise values to these bare parameters, since the 
observable parameters are those of the effective model 
to be derived below. They are the important ones, and 
they are best taken from experiment.) With t p d a ~ 1.80 
eV, only the antibonding band crosses the Fermi energy. 
At the Fermi energy, the wavefunctions have substantial 
weight at both copper and oxygen sites. 
With interaction, we must add a term 



Hi„t = U d ndr^ridrl 



-U p ^2n pr ^n P rii. (10) 



Here the sum runs over all atoms in the plane, and the 
interaction is approximated to be local; only electrons on 
the same atom interact. Theories of high-temperature 
superconductivity based on charge fluctuation-mediated 
attractions go beyond this approximation by includ- 
ing longer-range terms, particularly interactions between 
electrons on neighboring copper and oxygen atoms. 
These interactions are not so important for the magnetic 
mechanisms which are the subject of this paper. 

Since the parameter U d is large, (~ 10 eV) we may 
expect a considerable revision of the energy levels of the 
noninteracting problem to take place. In the ground state 
of the half-filled band case, the copper orbital will have 
essentially only one electron, since the addition of a sec- 
ond would cost this very large amount of energy. This 
leads to an insulating state because any flow of charge 
would require a real change of the occupations from their 
base values of one per copper and two per oxygen. Such 
a state is in fact observed in La2Cu04. It is also a mag- 
netic state, in the sense that every copper atom has spin 
one-half. (The atomic configuration is 3d 9 .) 

U p is smaller than U d , but large enough that holes 
added to the half-filled band state go predominantly to 
the oxygen sites. Each such site has then also a spin one- 
half. It was shown by Zhang and Rice |2§] that these 
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oxygen spins pair with the copper spins to form a spin 
singlet. This is the entity which moves from unit cell 
to unit cell. It produces metallic conduction, as it has 
a charge of —1 relative to the rest of the lattice. It also 
has a spin of zero, while the other cells have spin one-half 
because of their unpaired copper spins. The Hamiltonian 
for these particles is 



H Z R = ~t ^ C la c r+S,a ~ *' £ 

+f/^n rT n ri . 



(11) 



We refer the reader to Ref. 26 for the details of the 
derivation. Here c rt7 creates a hole and r refers to the po- 
sition of a unit cell, not an atom, and S and 7 are nearest- 
neighbor and next-nearest-neighbor vectors, respectively, 
of the square lattice of unit cells. This Hamiltonian is 
known as the single-band Hubbard Hamiltonian. It has 
received a great deal of attention in connection with the 
high-T c problem. Our aim here is not to review this work, 
but to do the simplest calculations which are relevant to 
superconductivity in a certain limit of the model, namely 
the limit U >> t,U >> t' . We shall assume that t and 
t' are of the same order of magnitude. 

The appropriate method in this limit is perturbative, 
in that the interaction term Hint is taken as the unper- 
turbed Hamiltonian. We have that 



Hi: 



and we wish to solve 



Ir-fTlrJ., 



Hint* - EM, 



(12) 



(13) 



to get started. Since H ln t simply counts the number of 
doubly-occupied sites, however, the basis of wavefunc- 
tions in configuration space (where the occupation of 
each site is specified) is already diagonal. The eigen- 
values are just 0, U, 2U, and so on, corresponding to 
0, 1, 2, etc., doubly-occupied sites. Each eigenvalue is 
very highly degenerate. Since our interest is in the low 
energy states, let us concentrate on the subspace with 
no doubly-occupied sites, with an unperturbed energy of 
zero. 

Now turn on the perturbation, which is 



(14) 



ks 



The effect of this operator on a state in configuration 
space is to move one electron from one site to a nearby 
site. In acting on a state with no doubly-occupied sites, 
it may, for example, move an electron to a site which is 
already singly-occupied, thereby increasing the number 
of doubly-occupied sites by one. The kinetic energy op- 
erator has the effect of mixing the different subspaces. 
We shall take this into account in the following way. If 



an arbitrary wavefunction is written in configuration 
space, then it may be decomposed as 



(15) 



where has no doubly-occupied sites, *f?i has one 
doubly-occupied site, and so on. 

The Schrodinger equation may then be written as 



/ Hqo 

Hio 



V 



01 





11 


H12 


21 


H22 









(*»\ 














= E 




) 


V : i 




I 1 / 



(16) 



Here the \1/ are as defined above, and H nn ' is the part 
of the Hamiltonian which acts only on a state with n' 
doubly-occupied sites and produces a state with exactly 
n doubly-occupied sites. This Schrodinger equation is 
still completely general. If U is large enough, and we are 
interested only in the low-lying states, we may neglect 
and all states with more than two doubly-occupied 
sites. Then the first two equations are 



(17) 
(18) 



Substituting the second equation into the first, we find 



Hoo*o + H 01 (E - H u )- l H 10 * = E^c 



(19) 



The second term on the left-hand side has the following 
structure. Wio^o is a state with exactly one doubly- 
occupied site. Its energy is therefore U, with corrections 
of order t or t'. H\\ acting on this state gives U to the 
approximation in which we are interested, and we may 
write the entire equation as 



Hoo^o - ^jHqiHw^o = EV . 



(20) 



We have then a new Schrodinger equation to solve, which 
may be written as 



with 



Hcs^a = E<$ 



Hcs — Hqo — —HqiHiq. 



(21) 



(22) 



In this equation, the second term represents virtual 
processes in which a doubly-occupied site is first created, 
then destroyed. In second order perturbation theory (or- 
der t 2 ), these processes affect the energies of the low- 
est energy states and break the degeneracy of the lowest 
level. Note that Hoq also has kinetic energy terms if the 
band is less than half filled, as it is possible then to move 
electrons from one singly-occupied site to another, cre- 
ating no doubly-occupied sites. Explicitly, the various 
operators have the form 
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Hio = -t 22 n r _ a ct !t7 c I+Sta (l - n r+s _ a ) 

—t n r, -cr C r, 0-^+7, cr(l — n r+7,-fr) (23) 

( X ~ n r,-J c L c r+ X( x"r+ 7 ,- CT (24) 
H °0 = -t (1 - ,- CT )4, ( r c r+5, -(l - "r+5,-<r) 

^ ( x ~ n r,- e r)4 )O -c r+0 ., o -(l ~ n r+0 , (25) 



The notation indicates that each nearest-neighbor pair 
and next-nearest-neighbor pair is counted once and once 
only in the sum. 

Writing the result out in its complete form, we have 



W e ff — t C r,er c r+<S,<7 t J]] c l,a 

r,8,a r,f,a 

J 

+ 2 



r.S 



r' L r+8 



J 



J J2 ( C l+8,a n r,-„< 



r+8' ,a 



C r+S,a C r,-cr C r+8' ,-o C r,cr 



(26) 



This rather complicated Hamiltonian can be separated 
into three physically quite distinct parts. The first term 
is a hoping term of the usual kind. An electron hops from 
the site r + 8 to the site r. The only point which must 
be borne in mind is that no doubly-occupied sites can 
be created in this process. This kinetic energy term is 
formally of order t. This is the leading term in an expan- 
sion in the small parameter t/U. However, the constraint 
of no doubly-occupied sites means that only "holes" (va- 
cant sites) can move. Thus, the contribution of this term 
to the energy is proportional to the doping level. If we 
call the density of vacant sites S, then the total contri- 
bution of this term to the energy is of order zSt, where 
z is the coordination number, (z = 4 for the square lat- 
tice). The next term is a spin-spin interaction. It lowers 
the energy of antifcrromagnetic spin configurations. The 
contribution of this term to the total energy is propor- 
tional to z(l - 8)J ~ 4z(l - S)t 2 /U. The final group of 
terms consists of the three-site terms, which represent a 
kind of induced kinetic energy through virtual processes. 
Their contribution to the energy is of order z 2 5J . In the 
high-T c materials, we normally have J/t w 1/3, and 6 
ranges from 0.0 to 0.4. As a function of 8, therefore, we 
expect the spin-spin interaction to be the dominant term 
at very small S, while at moderate values of 8, perhaps 
8 w 0.1, the other two terms will start to become impor- 
tant, based on this simple analysis. In fact, which terms 



are important depends a good deal on the property of 
the system in which we are interested. 



IV. THE DILUTE LIMIT 

The t-J model has strong support as a description of 
the low-lying energy states of the doped CuC>2 planes 
in the high-temperature superconductors. Such an im- 
portant model must be investigated over its full range 
of parameters, not only for the sake of understanding a 
model which is interesting in and of itself, but also for 
the purpose of gaining insight into the various trends of 
the model as its many variables change from one area 
of parameter space to another. In other words, such a 
wide-range investigation of the t-J model could link non- 
physical, yet soluble limits of the model to more physical, 
insoluble regions. In this spirit of identifying trends, we 
present in this section an analytical solution to the di- 
lute or low spin-density limit. The basic t-J model will 
first be considered, followed by the model with the ad- 
dition of three-site terms ]2?| and next-nearest-neighbor 
hopping. Justification for this "extended" t-J model is 
made below. 



A. The basic t-J model 

Consider two spins, one up and one down, on a square, 
two-dimensional lattice of N s sites. The goal is to find 
the ground state of the system satisfying the Schrodinger 
equation 



H tJ \iP)=E\^) 



(27) 



where E is the binding energy of the two spins. These 
spins interact as described by the "traditional" t-J model 
Eq. (Eq) which may be written as: 



H 



tJ 



(28) 



<j n r+8,a + ( 4+ < 5,(T C r,-CT C r+<5, -<r c r,cr) 



V,8,(T 



where summations are over the N s sites r, the four 
nearest-neighbor sites 8, and spin up and down a. The 
J-term is written in the above form for convenience in 
comparing it with additional terms later. To find the 
ground state of the system, the following pair wavefunc- 
tion form is used: 



riT C r 2 | 



10), 



(29) 



where the solution entails finding a(r) at all r. Before 
this is done, it should be noted that, although the t-J 
model implicitly includes the restriction of no-doubly- 
occupied sites, the given form in Eq. (pq) will not reflect 
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this constraint unless an on-site repulsion term is added 



T~ttj — T~(-tj + Hv , 



(30) 



where the limit Vq — > oo will yield the desired constraint 
a(r = 0) = 0. In other words, the probability of finding 
the pair on the same site will be zero as required. This 
added term is, of course, identical to the Hubbard [/-term 
of Eq. @. 

Substitution of Eqs. (§§) and © into Eq. @ and 
appropriate relabeling of summation variables reduce the 
problem to the following equation in a: 

- 2t^ a(r + S) - ~ 5>(r) - o(r - 25))S r . s 
s 1 s 

+V a(r)6 rfi = Ea{r), (31) 

where ri and r2 are fixed, and the relative coordinate 
r = ri — r*2 is used. Using the following relations: 



a(r) = ^ e lkr a(k), 



o(r), 



(32) 
(33) 



the Fourier transform of Eq. ( |3l| ) is taken, and solving 
for a(k) yields 



n \ 3 ST <x\ cosk ' 6 



e(k) - £7/2 
1 



'2N S Q{ j e(k)-S/2' 



(34) 



where e(k) = — 2t(cosfc x + cosk y ) is the tight-binding 
ener gy dispersion. The inverse Fourier transform of 
Eq. (B4) is then taken, resulting in the following: 



a(r) 



Jv^ fx s 1 e lk r cos k • 6 



s 

2 



N„ ^ e (k) - £/2 



V, ^e(k)-E/2' 



(35) 



a(0) can be found self-consistently from Eq. (35) to be 
J h 



o(0) 



v Jo + 2/y 



i 



iV, ^ e(k) 



£7/2' 



1 ^ — > cos h x 
X = W S *-' e{\l)-E/2- 



(36) 
(37) 
(38) 



Substituting the solution of a(0) back into Eq. ( p5| ) then 
yields the following: 



a{v)= J -Y,a{8) 
s 

1 v-^ e l 



(cosk-5-7 1 /(7 + 2/^o)) 



s k 



■ (k) - £7/2 



(39) 



It can be seen readily at this point that o(— r) = a(r) 
as should be expected from the symmetric form of the 
wavefunction chosen in Eq. (^9j) above. The four nearest- 
neighbor a(<5)'s can now be found. Letting 8\ = (1,0), 
^2 = (0, 1), £3 = (—1, 0), and 64 — (0, —1), the even sym- 
metry of a(r) requires that 0(^3) = a(Si) and 0(^4) = 
0(^2)- Solving for both a(Si) and 0,(62) self-consistently 
in Eq. ( |39| ) results in 0(62) — ±a(5i). Thus, two solu- 
tions for a(r) are obtained: 

1 



s k 



(cosfc^ + cosky - 2I 1 /(I + 2/V )) 



e(k) - £7/2 



a d (r) = Ja(<5] 



N s ^ 

s k 



. (cos k x — cos k y 
e(k) - £7/2 



(40) 
(41) 



where Eqs. ( |40| ) and (fll]) correspond to an extended s- 
wave (s ;E 2 + j / 2) solution and a d-wave (d x 2_ y 2) solution, 
respectively. 

The (4,2 _j,2 -wave solution in Eq. ( [i"l| ) has no depen- 
dence on Vo, and it can be seen by symmetry alone that 
a d(0) = 0; that is, there is no double-occupancy. The 
s-wave solution in Eq. ( fh]| ) also obtains this constraint 
upon letting the on-site repulsion Vq become infinite: 



Vn 



lim a s (r) 



s k 

(cos k x + cos k y — 2 Ji / Jq ) 
X e(k) - £7/2 



(42) 



where it can be seen that a JO) = by using the defini- 
tions of Iq and I\ in Eqs. ( |37| ) and (38) above. Finally, 
in the thermodynamic limit (N s — > 00), the sums over k 
become integrals over <ik, and the solutions to the "tra- 
ditional" t-J model in the dilute limit of one pair of spins 
become 



a s (r) = Ja(6 x ) 



1 



dke 1 



(2irf 

(cos fca; + cos k y — 2Ix/I ) 



a d (r) = Ja(5i) 



e(k) - £7/2 



dke 1 



/o = 



/1 = 



(2tt) 2 

(cos k x — cos fcj,) 
e(k) - £7/2 

1 'dk 1 
dk 



(27T) 2 
1 

(2^ 



£ (k)-£7/2' 

cos k x 
e(k)-£7/2' 



(43) 

(44) 
(45) 
(46) 
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where integration is over the full Brillouin zone: k x = —ix 
to 7r and k y = —it to w. 

Of some interest is the critical value of J at which 
there is a bound state for each solution. Letting r = <5i 
in Eqs. ([l3]) and (^3), the relationships between J and E 
for the two symmetries can be obtained: 



1 



1 



(2tt)2 

1 



^ cos k x (cos k x + cos k y - 2h/I Q ) 



Jd (2tt)* 



dk 



e(k) - E/2 
cos k x (cos fcj; — cos k y ) 
e(k) - £/2 



(48) 



These integrals cannot be done in closed form but 
lend themselves fairly easily to straightforward numeri- 
cal methods such as Monte Carlo integration. Some care 
must be taken when considering small E where the nu- 
merator tends to zero at k = 0. In the J — ► oo limit, 
there are precisely four bound states corresponding to 
the four sites on which the attraction is nonzero. Two of 
these are p- wave, and therefore not allowable. Thus there 
are only two bound state solutions in this limit. Since the 
number of bound states can only decrease as the attrac- 
tion is decreased, the maximum number of bound states 
is two. The s-wave and d-wave solutions discussed here 
are the only possible bound states in the model. 

Equations (47) and (E3) are the basic equations deter- 



mining the bound state energies, and we derived them 
from the Schrodinger equation. However, we may derive 
them from the BCS gap equation as well. In this re- 
gard, see also Refs. p^-|30[|. At zero temperature the gap 
equation is 



A(k) = -£V(k,k') 



A(k') 



vwH*F+A(kp 



(49) 



In the limit when |A(k)| is small, this becomes linear: 

A(k) = -X> ( k,k0^2_ (50) 



This is the stability form of the gap equation. If it has 
a solution, then the system is superconducting (with in- 
finitesimal gap) at zero temperature. The reason this 
form of the gap equation is rarely seen is that, in conven- 
tional superconductivity with BCS-like purely attractive 
interactions, superconductivity at zero temperature oc- 
curs even with infinitesimal interaction strength. Here 
this is not necessarily the case. Note also that no fre- 
quency cutoff is assumed in the interaction, because the 
basic model contains no retardation. This is another im- 
portant difference from the usual BCS model. 

We may now transform Eq. ( |50|) into Eq. ( [47| ) by 
first making the substitutions (i — 0, (corresponding 
to the dilute case), A(k) = A^cos^) + cos(fc y )), and 
V(k, k') = — J(cos(k x — k' x ) + cos(k y — k' )), then multi- 
plying the resulting equation by cos(k x ) + cos(k y ) and in- 
tegrating over k. The result is Eq. (Effl . A similar proce- 
dure using A(k) = Ad(cos(A; x ) — cos[k y )) yields Eq. (48). 
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Jit 

FIG. 2. Phase diagram of binding energy E vs. coupling 
constant J for the extended t-J model in the dilute limit 
for (a) t' = 0.00 t, (b) t' = -0.15 t, (c) t' = -0.30 t and 
(d) t' = —0.45 t. The solid line plots correspond to calcula- 
tions done without the three-site term (J3 = 0), the dashed 
lines for those done with the three-site term (J3 = J), the 
plots with open diamonds are for s ;c 2 +!/ 2-wave, and the plots 
with closed squares are for d a .2_ H 2-wave. The addition of the 
J3-term causes the s-wave curves to shift left, whereas the 
d-wave curves shift right off to infinity; that is, there is no 
d-wave solution. The addition of the t'-term causes d-wave 
(for J3 = 0) to shift left and cross the s-wave curve at a point 
dependent upon the strength of t' . 
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Hence the gap equation at zero gap is the same as the 
Schrodinger equation at zero binding energy. 

A phase diagram of E vs. J depicting the s- and d- 
wave bound states may be seen in Fig. 0(a). It can be 
noted from this low density phase diagram that s-wave 
always has the lower energy of the two solutions at fixed 
J and is, therefore, the ground state of the problem at 
hand. Furthermore, s-wave has a bound state above a 
critical- J ( J c ) of 2.00 t whereas d-wave has a J c of about 
7.32 t. The phase diagram of Fig. ||(a) is identical to that 
reported recently by Hellberg and Manousakis j3l[] and 
the s-wave results are also similar to those of Kagan and 
Rice @. 

If one holds to the idea that magnetic pairing should 
produce d-wave, because of the hard core in the poten- 
tial, the surprising aspect of these results is that s-wave 
always has by far the largest binding energy of the two 
symmetries. We may in fact prove a general theorem that 
if the potential has square symmetry, the ground state is 
always s-wave for nearest-neighbor hopping. This theo- 
rem holds whether there is a hard core or not. The proof 
is given first on the continuum with a radially symmetric 
potential, since it illustrates the main ideas most clearly. 
The Hamiltonian is 



H 



2(i 



ld_ 

pdp 



d 
% 



1 d 2 

72 PiAZ 



+ V{ P ). (51) 



We argue by contradiction. Let ip m (p, 4>) be the normal- 
ized ground state and have non-s-wave symmetry: 



1 



f(p)e 



im<f> 



(52) 



with to ^ and J Q pf 2 {p)dp = 1. The ground state 
energy is 



(53) 



2p J dp \ dp 

m 2 h 2 r 00 i , 2 , . , 

■ / -J (P)dp 
Jo P 

> 

pf{p)V{p)dp. 



2[i 



Now consider the normalized trial wavefunction 
ip s (p,4>) = f(p), which is s-wave. This has the expec- 
tation value 

E s = (il> s \H\ip s ) 



pf 2 ( P )V(p)dp. 



Clearly, 



E s < Ed — E s + 



2p Jo P 



: f(p)dp. 



(54) 



(55) 



Thus ip m cannot be the ground state. A similar proof ob- 
viously holds in three dimensions. The physics is simply 
that the increased kinetic energy from angular variations 
in a higher- wave wavefunction will always be greater than 
any gain in potential energy coming from avoidance of the 
core. 

On the lattice, the physics is clearly the same, and 
the proof is quite similar. Let the potential V(r) have 
square symmetry, and let i r . r ', the hopping coefficients, 
be positive. Let the assumed ground state be tpd( r ) have 
d-wave (-Bi) symmetry. It may be taken to be real. Now 
consider the wavefunction 



V's(r) = \ipd( r )\- 
The expectation values are: 

E d = (Mn\ip d ) 



t r ,v>ipd(r)ipd(r') 



(56) 



(57) 



and 



E s = (il> a \H\il>.) 

= ]Tv»^(if 



tr^\Mr)M^)\- (58) 



Comparison shows that 

E s < Ed, 



(59) 



so that tpd can never be the nondegenerate ground state. 
Except for pathological cases, it is not hard to improve 
ip a to get a strict inequality in Eq. [59]. (Let ip s (x,x) — > 
i/j s (x, x) + 6i/j, with 5 sufficiently small). We have there- 
fore shown that s-wave is in fact always the lowest energy 
solution. Hard core or not, <i-wave never wins. 



B. The addition of the three-site terms 

We now include the so-called three-site terms of the 
t-J model in our dilute limit investigations. These added 
terms originate from the derivation of the t-J model from 
the Hubbard model. Considering the Hubbard model, 



bbard 



(60) 



the t-term is treated as a perturbation with respect to the 
J7-term, where U t. Only terms up to and including 
order t 2 /U are retained. Letting J = At 2 /U , the pertur- 
bative expansion yields the t-J Hamiltonian of Eq. ( psj ) 
plus the three site terms |27j . The Hamiltonian may now 
be written 
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H 



uj 3 



(61) 



J: 



<5,ct "T" C r+<5,<T C r,-o- c r+<5,-<T C r,o-) 



r, 5' ^S.a 



+ C l+5,A-<T C r+5>,-* C r,*)i ( 62 ) 

where J3 is merely used as an adjustable factor; letting 
J3 = restores Eq. (^8|), whereas letting J 3 = J yields 
the extended model as derived from the Hubbard model. 

The three-site terms are usually dropped due to their 
complexity and the fact that their expectation value is 
formally of order 5 J, not J, where 6 is the hole den- 
sity. However, they are enhanced by the fact they are 
of order z 2 , not z, where z is the coordination number. 
Numerical results show clearly that they become impor- 
tant when S is greater than 0.05 Certainly, it would 
then be expected in the dilute limit that the three-site 
terms become very important interactions concerning the 
ground state energy of the system, and the solution to 
the extended model in the two-spin limit bears this out. 

The solution of two spins interacting on a two-dimen- 
sion al squ are lattice of N s sites is identical to that of 
Sec. IV A and need not be repeated here as now applied to 
Eq. (|6l|). The wavefunction coefficients a(r) of Eq. ( |29| ) 
are found to have the following solutions (upon letting 
N s -» 00) 



a s (r) = j(l + ^)a(J 1 )^ 

(cos + cos ky — 2Ji/J ) 



dke lkr 



«d(r) 



e(k) - E/2 
(cos k x — cos ky) 



dke 



ik-r 



■ (k) - E/2 



(63) 



(64) 



where Io and I\ are still defined as in Eqs. fl45|) and ( |46[ ), 
respectively. As expected, when J3 = 0, the solutions 
of the t-J model without three-site terms, Eqs. (|^|) and 
(p4|), are retrieved. 

Given that J3 = J, it can be seen that there is only 
a trivial d-wave solution, namely cid(r) = for all r. In 
other words, there is no bound state solution for d-wave 



pairing. The s-wave solution is similar to that of Eq. ( |43[) 
differing only in pre-factors. Finding the relationship be- 
tween J and E for these two solutions yields 



1 



dk cos k r 



J J (2tt) 2 
cos k x + cos k y — 2Ii I To) 



Jd 



e(k) - £/2 
J 3 \ 1 



(65) 



J / (2tt) 2 



dk cos /c T 



(COS fca, — COS k y ) 

e(k) - E/2 ' 



(66) 



It can be readily seen that, as the three-site interactions 
are "turned on" , that is as J3 is increased from to J, 
the J vs. E curve for <i-wave would move off to infinity 
in the phase diagram of Fig. ^(a). This again shows that 
there is no bound state solution for d-wave in the dilute 
limit of the extended t-J model. For s-wave, the J vs. E 
curve is exactly one fourth that of the s-wave solution 
for the "traditional" t-J model as shown in Fig. ||(a). 
Thus, the three-site terms are s-wave enhancing, <i-wave 
suppressing interactions. 

The expectation values of the various terms of the ex- 
tended t-J model can be expressed in terms of the wave- 
function coefficients a(r) as well as in integral forms. The 
normalization is 



(v#) = i a ( ri ~ r2 ) 



(2-k? 



dkAl 



(67) 



where N s is temporarily left finite but will cancel out in 
the expectation values shown below. The coefficients Bj a 
and Ak are defined for convenience as 

A k = (cosfc, + cosky - 2/ 1 // )/(£(k) - E/2), 
B, h = 1 + 3J 3 /J (68) 

for s-wave and 

A k = (cosk x - cosfc. y )/(e(k) - E/2), 
B Js =l- J 3 /J (69) 

for d-wave. 

The expectation values can now be expressed in terms 
of integrals. The relations between J and E from 
Eqs. (|6"5|) and (|66]) are also employed to yield 



(Ht 



(</#> 

J dkA k (cosfc x + cosky) 



-4t 



JdkAl 



<Hj a ) = 



4 J dkAk cos k x 
JdkAl ' 

12_ /J3 

Bj 3 { J 

+^ 
B.; 3 \ J 



J dkAk cos k x 
JdkAl ' 
J dkAk cos k x 
JdkAl ' 



(70) 
(71) 

s-wave (72) 
d-wave. (73) 



Of course, the solutions for d-wave are ill-defined for J 3 = 
J (dd(k) = for all k in this case) but are included here 
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merely for completeness. What can be immediately seen 
is that, for s-wave, the three-site term has an expectation 
value exactly three times that of the J-term, or (Tij 3 ) = 
3(7ij), and for <i-wave, the two expectation values are 
equal and opposite, or (Ti.j 3 ) = — (Hj)- 



/-term interactions 
Initial State Virtual State 




Final States 



/ 3 -term interactions 
(three possible processes) 



Initial State 



Virtual State 




FIG. 3. Diagram of the dynamics of the spin-spin corre- 
lation term (or J-term) and the three-site term (or J3-term). 
Shown are initial, intermediate (or virtual), and final states 
for one example of each interaction (filled arrows). For the 
J3-term, alternative interactions are also shown (open ar- 
rows). In these examples, a down spin has a virtual hop- 
ping to a nearest-neighbor site occupied by an up spin. In 
the J-term, either spin hops back to the original site, and in 
the Js-term, either spin hops to one of the other three near- 
est-neighbor sites. In the dilute limit, the spins have equal 
probability of hopping to any of the four nearest-neighbor 
sites as shown here. In the Cooper problem, however, there 
would be a probability of 5 for each of the three originally 
unoccupied nearest-neighbor sites to actually be occupied by 
a non-interacting spin. 



A closer look at the various interactions of the ex- 
tended t-J model in the dilute limit is worth making at 
this point. It is obvious from the Hamiltonian as written 
in Eq. (^lj) that the J3-term is quite similar to the J- 
term. The three-site terms comprise interactions involv- 
ing a site and two distinct nearest-neighbors, whereas the 
J-term is the special case where the two nearest-neighbor 
sites are merely the same site. Figure || shows a depiction 
of the virtual processes involved in the spin-spin corre- 
lation term (J-term) and the three-site terms (J3-term). 
When the two spins are on nearest-neighbor sites, the 
J-term can be thought of as a virtual hopping of one 
spin onto the site occupied by the other spin, followed by 



either spin hopping back to the first site. In the dilute 
limit, the probability of finding the other three nearest- 
neighbor sites unoccupied is exactly 1. Hence, when one 
spin has a virtual hopping onto the site occupied by the 
other spin, there is equal chance of one of the spins hop- 
ping onto any of the four nearest-neighbors. For s-wave, 
it is easy to see then why the J3-term is three times the 
size of the J-term, and, due to its change of sign with tt/2 
rotation, the d-wave Jj-term is equal to the J-term but 
with opposite sign. Hence, the total summation of the 
J3- and J-tcrms for <i-wave yields zero due to its symme- 
try, and the total summation for s-wave yields four times 
the value of the J-term. The importance of this result 
will become more apparent in Sec. ^ where increasing 
spin density is considered in the Cooper problem. The 
actual effect of the three-site terms is more akin to a 
kinetic energy term than an interaction term. It is an 
effective hopping to second- and third-nearest-neighbors. 
It is therefore not so surprising that it tends to stabilize 
the s-wave state. 



C. The addition of the next-nearest-neighbor 
hopping term 

One more interaction is included in our extended t-J 
model, namely the next-nearest-neighbor hopping term. 
This term is added with the intent of making a first ap- 
proximation at modeling specific high-temperature su- 
perconducting systems. It is known, for example from 
neutron scattering, that the various high-T c copper- 
oxides differ in energy dispersion, Fermi surface, etc. Us- 
ing the next-nearest-neighbor hopping term with coeffi- 
cient t' is a way of distinguishing these compounds from 
one another in our tight-binding model. Perhaps a more 
systematic approach, emulating our derivation of the J3- 
term, would be to begin with the Hubbard model and 
include the i'-term in the large-t/ perturbative expan- 
sion. This would result in more terms of order J, that is 
various three-site terms involving nearest-neighbor sites, 
next-nearest-neighbor sites, and mixtures of both. The 
many extra terms may be worth studying at some future 
time, but for the purposes of identifying trends in the t-J 
model, we do not wish to introduce even more complex- 
ity to our extended model until it has been sufficiently 
investigated in its present form. 

The extended t-J model is given in its full form as 



Utt'JJs ~ 



t ^ t c r,er c r+<5,<r 
J 



* E 4,. 



(74) 



■<r' t r+<5,<7 



J3 

4 



+ C r+S,a C l,-a c r+S,-a c r,a) 
( C l+S,cr n r,-cr C r+ S' .cr 
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i-6> C r, 



ho ,— <7 



where summation over 7 means summation over next- 
nearest-neighbor sites. Once again considering the prob- 
lem in the dilute limit, the solution to the problem of 
two spins interacting with the Hamiltonian of Eq. (|74|) 
results in exactly the same forms for a s (r) and a<j(r) 
as given in Eqs. (|6^) and (Q), respectively, only now, 
the tight-binding energy dispersion is given by e(k) = 
— 2i(cos k x + cos k y ) — At' cos k x cos k y . 

The relationships between J and E are of the same 
forms as given in Eqs. (|65| ) and (|66| ) for s-wave and d- 
wave, respectively. Once again, phase diagrams for E 
vs. J can be drawn for various values of t' . Fig.|2j(a) 
is, obviously, the phase diagram with t' = 0. Figs. 0(b), 
(c), and (d) show the phase diagrams for t' = —0.15 t, 
—0.30 t, and —0.45 t, respectively. It can be seen that, 
for finite t' and no three-site interactions, the s-wave and 
e?-wave curves cross at a binding energy E and value 
J which decrease with increasing magnitude of t' . In 
other words, as the next-nearest-neighbor hopping term 
increases in strength, there is a cross-over between an 
s-wave ground state solution and a d-wave ground state 
solution at smaller and smaller binding energies and J- 
values. This shows that the i'-term with a negative sign 
is an s-wave suppressing, c£-wave enhancing term in the 
model. That this must be the case is evident already 
from Eqs. (|57j) and (^8|). Once again, when the three-site 
terms are considered (J 3 = J), there is no d-wave solu- 
tion in the dilute limit, and the s-wave curve is reduced 
to a quarter of its value without the J3-term. 

The physics of the dilute limit is that, in the absence 
of a Fermi surface, the s-wave instability is dominant. 
The s-wave simply has lower kinetic energy. In the next 
section we present an investigation of a single pair now 
interacting in the presence of a Fermi surface: the ex- 
tended t-J model is evaluated in the Cooper problem. 



V. THE COOPER PROBLEM 



The investigations of Sec. [V can now be extended from 



the dilute limit to the regime of finite density. Again, the 
purpose of this study is to identify trends with the ad- 
dition of the three-site term, with the changing of the 
next-nearest-neighbor hopping term, and now with the 
changing of the density of spins. As before, the prob- 
lem is solved for one pair of interacting spins, now being 
in the presence of a filled Fermi sea of non-interacting 
spins. These results will be compared qualitatively with 
previous VMC results which also yield some of the trends 
detailed in this work. 

Consider the problem again of two spins interacting 
on a square, two-dimensional lattice of N s sites, but now 
allow them to interact in the presence of a filled Fermi 
sea up to a Fermi energy (chemical potential) of /x. This 
is the Cooper problem as applied to the extended t-J 
model. As we shall show, this also corresponds again to 



the stability form of the gap equation at zero temperature 
with a certain interaction. 

The solution once again entails finding the wavefunc- 
tion coefficients a(r) Eq. (^9|) from the Schrodinger equa- 
tion Eq. (p^) which is now rewritten 



H 



tt'jj 3 



(2a* + E) 



(75) 



where the eigenvalue E has been redefined with respect 
to the Fermi energy. Once again, Htt'jj a must include 
an on-site repulsion term so that the solution will enforce 
no-double-occupancy as was shown in Eq. ( |3C| ) . In a sim- 
ilar fashion, the states within the Fermi surface must be 
excluded from any interactions since they are to be filled 
with non-interacting spins. This is still a two-particle 
problem which reduces to a single-particle equation as in 
Eq. (|3l|), but with a fictitious interaction coming from 
the Pauli exclusion principle. The exclusion of states 
within the Fermi surface is taken into account by adding 
another term to Ti.tt'jj 3 , namely 



T~ttt'jj 3 

with the following dependence on k 



k,o- 



(76) 



where £(k) = e(k) — ji. Thus, the same technique as with 
the exclusion of doubly-occupied sites can be employed: 
the limit of Uq — ► 00 at the end of the derivation will 
exclude states within the Fermi surface. 

An important observation needs to be made concern- 
ing the three-site terms. Figure shows the dynamics 
of the J-term and J3-term in the dilute limit, where the 
J3-term has three times the value of the J-term for ex- 
tended s-wave and the same value and opposite sign for 
d-wave. This is due to the fact that, if the two spins 
occupy nearest-neighbor sites, the other three nearest- 
neighbor sites would then be guaranteed unoccupied in 
the dilute limit. Now, however, with finite density, if the 
two interacting spins are on nearest-neighbor sites, then 
there would be a virtual hopping onto one of the sites, 
and the probability of the original nearest-neighbor being 
unoccupied would be 1, whereas the probability of the 
other three nearest-neighbors being unoccupied would 
be 8, the hole density. In other words, the three-site 
terms need to have a pre-factor of 5 to take into account 
the finite density of non-interacting spins now present on 
the lattice. In addition to this, the nearest-neighbor and 
next-nearest-neighbor hopping terms must also both have 
pre-factors of 5 to reflect this finite probability that a site 
to which one of the two interacting spins may hop could 
be occupied by one of the non-interacting spins. Given 
all of these necessary changes, the Hamiltonian (without 
the additional Vq- and £/c-potential terms) should now 
be expressed as 
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(78) 



"4 2 ( n r,- CT «r+5. C 
r,<5,cr 



-4 E 



CT n r,-cr C r+<5' ,<r 



,o- C i\-<x C r+<5',-er C r,er)j 



Upon solving for a(r) and self-consistently solving for 
a(0) and the four a(5)'s, two solutions are again obtained 
of the forms 



a s (r) = j(l + 3S^\ ^1)7^2 / dkelk r 

(cosfe^ + cos ky — 2I\/Iq) 
X S£(k)-E/2 + U c (k) ' 

a„(r) = J (l - 4) a (*i) (2^)2 / dke *" 

(cosk x — COS fcy) 
X <S£(k)-i;/2 + [/c(k)' 



(79) 



(80) 



which once again correspond to extended s-wave and d- 
wave solutions, respectively.. Here, the limit Vq — > 00 
has already been done. Now the definition of Uc(k) from 
Eq. ( fn| ) can be used, and the limit of Uq — > 00 is taken, 
yielding the final solutions 



oa(r) = J ( 1 + 3^40 a (<5i)7^ / dke*" 



(cos fe x + cos fe a — 2Ji /Jo) 
J£\ . 1 

(cos fcj, — COS fey) 
<^(k)-£/2 ■ 



(81) 



o«,(r) = J(l-*-i)a(5i) 



ik-r 



C(k)>0 



(82) 



Integrations are, therefore, done over those states lying 
outside of the Fermi surface; the filled states do not con- 
tribute to the interactions of the problem. 

The relationship for J and E can be found for the 
finite density problem as well as the expectation values 
of the various terms of the extended t-J Hamiltonian. 
The equations connecting J and E can be expressed as 



1 

J 

where now 



B, 



1 



(2nf 



dk cos k x A k , 



(83) 



?(k)>0 



A k = (cos k x + cos k v - 2/ 1 /I )/(<J£(k) - E/2), 
Bj 3 = 1 + 35 J 3 / J (84) 

for s-wave and 



A u = (cos k x - cos k y )/(S^k) - E/2), 
Bj 3 = 1 - SJ 3 /J 



(85) 



for d-wave. The integrals in Eq. (|8^) cannot be done ana- 
lytically. The actual calculations entailed a simple Monte 
Carlo integration method where special care was taken 
when E was small and the denominator of the integrand 
tended towards zero along the Fermi surface. 

These eigenvalue-type equations can again be shown 
to be equivalent to the stability form of the gap equation 
at zero temperature. The difference with the dilute case 
is that the interaction V^ a p which appears in the gap 
equation now has a frequency cutoff 



V fk k"i-J y ( k ' k ')> £( k ) >0 
*Wk,k J - I Qj f(k)<0. 



(86) 



This looks artificial. However, the Cooper bound state 
equation in the conventional case has the same asym- 
metrical cutoff. The difference between the asymmetri- 
cal cutoff and the more symmetrical cutoff of the BCS 
interaction is in fact responsible for the factor of two dif- 
ference in the exponential factors in the Cooper binding 
energy and the BCS gap. Again, this difference does not 
influence the stability issue of interest here. 
The normalization is found to be 

= K r i- r 2)l 2 

ri,r 2 

= J^Bl\a(S 1 )\ 2 7 ^ [ dkAt, (87) 
and the expectation values of the terms may be expressed 



(n t ) = 



(V#> 

Je(k)>o d^ A l(cos k x + cos k y ) 



-ASt 



X;(k)>O dk ^k 



= -8St 



Jf(k)>o dk ^k 



(88) 



(89) 



(Hj) = 



(K, 



4 J m>0 dkA k cosk x 
12 f,J 3 \ J m>0 dkA k cosk x 



B Js V J J I m >o dkA l 



(90) 



,s, (91) 



4 / ,h\ Liu^n dkA k cos k x 



B -h V J J J iW >o dkA 
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FIG. 4. Phase diagram of coupling constant J vs. 
spin density n in the Cooper problem for binding energy 
E — —0.01 t, J3 = (solid line), and J3 = J (dashed line). 
The cross-over points between d- (closed squares) and s-wave 
(open diamonds) occur at specific densities for both cases, 
where these densities remain fairly fixed as E (and accord- 
ingly J) is decreased even more. With the addition of the 
three-site term, s-wave commands more of the phase diagram 
while d-wave is suppressed. 



Equation (|83j) is divergent in the limit E —> 0. In other 
words, given the presence of a Fermi surface (S < 1 and 
n > 0), J does not have a finite critical value as it did 
in the dilute limit (5 — 1 and n = 0). This means that 
there exists a bound state (E < 0) for any potential J 
no matter how small. This is indeed the conclusion of 
the Cooper problem: any pairing potential has a bound 
state for one pair of spins, and this must extend to the 
case of multiple pairs of spins. It is now important to 
see which of the two symmetries has the largest binding 
energy E (in absolute value) for a given J as both t' 
and 6 are varied. In reality, J is expressed as an integral 
containing E, and so E is fixed as J is calculated for 
varying values of both t' and 5. Fig. ^ shows the diagram 
of J vs. n for a fixed binding energy E — —0.01 t for 
t' = and J3 = 0. The curves for s-wave and d- wave 
cross at a density of about n — 0.2. This means that 
d-wave has the larger binding energy at smaller potential 
J for densities n > 0.2, whereas s-wave has the larger 
binding energy at smaller J for n < 0.2. Turning this 
statement around, given a set J, d-wave would have the 
higher binding energy for n > 0.2 and s-wave would be 
dominant for n < 0.2. Indeed, if the curves of J vs. n (as 
calculated from the above derivations) were to be plotted 
as in Fig. || for decreasing values of the binding energy 
(E — > 0), this cross-over point would remain relatively 
fixed even as the J-curves collapsed to J = over all 
n at E = 0. This cross-over between an s- and d-wave 
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FIG. 5. Phase diagram of next-nearest-neighbor hopping 
parameter t' vs. spin density n in the Cooper problem for (a) 
J3 = and (b) J3 = J. The phase boundary represents small 
binding energy E (and, therefore, small coupling constant J). 
These diagrams show that the f'-term enhances d-wave over 
s-wave, whereas the J3-term enhances s-wave over d-wave. 



phase near n = 0.2 has been seen by Dagotto, et ai, in 
a quite different low-density study ]33| . 

The interesting trends of these calculations can be seen 
when the various parameters are changed. When the 
three-site term is added, or J3 = J, this cross-over point 
shifts to the right in the phase diagram as seen in Fig. |] 
to around n — 0.53. Here, s-wave commands more of the 
phase diagram while d-wave is suppressed. This is con- 



sistent with the previous results of Sec. [TV] where the 
three-site term was found to be an s-wave enhancing 
interaction. The cross-over points were also found for 
various values of t'. These results are shown in Fig. ||. 
Here, the cross-over from d-wave to s-wave with decreas- 
ing n (increasing 6) occurs at lower n as the £'-term in- 
creases in strength. Once again it is seen that the next- 
nearest-neighbor hopping term enhances d-wave, where 
this symmetry takes over more of the phase diagram as 
t' increases in magnitude. It can also be seen even more 
clearly in Fig. [5] that the addition of the three-site term 
increases the size of the s-wave phase as has now been 
described in great detail. 
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The solution of the Cooper problem shows that the 
presence of a Fermi surface can stabilize d-wave over s- 
wave. To see this in more detail, note that the interaction 
in momentum space is V(k, k') = — (3/4)(cos(fc x — k' x ) + 
cos(fcj, — k' )) and the potential energy is proportional to 



]Tv(k,k>(k)a(k') 



(93) 



k,k' 



Now divide momentum-pair space into three regions: re- 
gion (I), k - k' « 0; region (II), k - k' m (±n,±n); 
region (III), k - k' « (2tt,0) or (0,2vr). In region (I), 
V(k,k') > 0, in region (II), V(k,k') < 0, and in re- 
gion (III), V(k, k') > 0. Furthermore in region (I), 
ad(k)a,d(k') > 0, in region (II), a^k^^k') < 0, and in 
region (III), eid(k)ad(k') > 0, for the d-wave case. For the 
s-wave case, a s (k)a s (k') > in all regions. In general, 
when regions (I) and (III) are important, then s-wave 
will have lower potential energy. When region (II) is im- 
portant, then d-wave has lower potential energy. The 
presence of the Fermi surface enhances the stability of 
the d-wave solution because it reduces the importance of 
region (I) and increases the importance of region (II). 

Although the dilute limit and the Cooper problem 
are not truly realistic pictures of superconductivity in 
the high-T c cuprates, these rather simplistic pictures do 
yield important behavior of the t-J model which extends 
into the many-body problem of multiple pairs of inter- 
acting spins. This crossing-over from d-wave to s-wave 
suggests regions of parameter space where d-wave and 
s-wave could even compete and form a kind of mixed s- 
d state. Section VI includes analysis of the mean field 



t-J model, which is shown indeed to support mixed s-d- 
pairing. 



VI. THE MEAN FIELD APPROACH 

A next logical step in the analysis of the extended t-J 
model is to evaluate its mean field Hamiltonian. As with 
the dilute limit and the Cooper problem, the goal is to 
gain physical insight into the contributions of the various 
terms of the model, but now by using a more complex 
picture than just a simple single-pair model. Of most 
interest is the form and symmetry of the gap function 
A(k) as obtained from the zero-temperature BCS gap 
equation. To derive A(k), the full Hamiltonian Eq. fl74| ) 
must first be reduced to the form 

H m { = ^2 e(k)nk, CT 

k,cr 

+ ^y(k,k04 T ct k; c_ ka c k , T . (94) 

k,k' 

(Throughout this discussion, J3 = J.) As with the 
Cooper problem, factors of 5 are included to account for 
the restriction of no-double-occupancy since Ti. m { does 
not take this explicitly into account. Equation (194), as 



applied to the mean field t-J model, is thus produced by 
transforming Eq. ( |T8| ) into a Hamiltonian of creation and 
annihilation operators in momentum space. e(k) is once 
again the tight-binding energy dispersion, and V(k, k') 
is of the form 

V(k, k') = -J(cos(k x ) cos(*4) + cos(k v ) cos(fc^)) (95) 
— 5J 3 [(cos(fc x ) cos(fc^) + cos(ky) cos(k' y )) 
+ 2 (cos(k x ) cos(k' y ) + cos(ky) cos(k x ))] . 

With V(k, k') now defined, H m f can be evaluated using 
the BCS wavef unction 



= Ufa + ^4 T c -ki) l°) > 



(96) 



where the BCS variational parameters are of the usual 
forms 



ul = li 1 - 



vl = r 1 



^ 2 (k) + A2(k) < 
Ve 2 (k) + A2(k) 



(97) 
(98) 



Here, £(k) = e(k) — /i and A(k) must be obtained from 
the zero-temperature gap equation Eq. (49). 

The gap function A(k) was calculated from the self- 
consistent BCS gap equation using an iterative numerical 
method. A matrix of initial gap function values was de- 
fined on a grid representing points in momentum space. 
A first iteration matrix was then calculated from the gap 
equation using the initial matrix as input. A linear com- 
bination of the first iteration and initial matrices was 
then defined and optimized with regards to the gap equa- 
tion, and this process was repeated until the matrices 
converged upon a steady solution. The integrity of the 
gap function solutions was checked (i) by using a variety 
of initial matrices, including matrices of random-valued 
elements, and (m) by changing the size of the matrices. 
The iterative method did indeed converge on reproduce- 
able, self-consistent solutions by using lattice sizes of 442 
to 1682 fc-states. 

The gap function for the mean field t-J model was 
found to always have three solutions which satisfied the 
self-consistent zero-temperature BCS gap equation. The 
first was the trivial solution A(k) = 0, and the other two 
were of the familiar forms 



A s (k) = ip s (cosk x + cosfc y ), 
A d (k) = ip d (cosk x - cos k y ), 



(99) 
(100) 



where the magnitudes ip s and tpd vary as functions of J, 
8 and t'. Under certain circumstances, there were also 
found to be s+d solutions when ijj s and tpd were com- 
parable in size, when J was larger than some threshold 
value of approximately 1.0 t, and when S was of some 
intermediate value 0.15 to 0.45. Hence, in at least a qual- 
itative fashion, sd-mixing is supported by the mean field 
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FIG. 6. Magnitudes of gap equation solutions to the BCS 
gap equation vs hole density for J3 = J = 0.33 t. Shown are 
(a) s x 2 +y 2- and (b) d a .2_ H 2-wave gap amplitudes for t' values 
of 0.00 t (solid lines), -0.15 t (dashed lines), -0.30 t (dotted 
lines), and —0.45 t (dot-dashed lines). The increase in ip s in 
(a) with t' is due to the flattening of the energy dispersion at 
low spin density (large hole density) as well as the narrowing 
of the band width as negative t' increases in value. The ipd 
peaks in (b) lie at hole dopings which closely match optimal 
doping in the materials modeled by the corresponding peaks' 
t' values. 



t-J model in regions of parameter space where the s- and 
d-wave phases are competitive. 

More interesting results were obtained from the mean 
field gap calculations where physical values of the varia- 
tional parameters were used. For J = 0.33 t, there were 
no significant mixed state solutions found for any values 
of 6 or t' , but the s- and d-wave solutions yielded sur- 
prising trends. Figure |^ shows the magnitudes of these 
solutions to the mean field gap equation over many hole 
densities and values of t'. The s-wave magnitudes are 
largest around lower spin density as expected from the 
Cooper problem results, however the values increase with 
more negative if . This trend of t' enhancing s-wave seems 
contrary to the previous results from the dilute limit and 
Cooper problem, but this incongruent result can be ex- 
plained by the nature of the energy dispersion. At these 
low spin densities, the filled states lie close the the en- 
ergy minimum. As t' becomes more negative, the energy 



dispersion becomes more fiat around the minimum, and 
the band gap reduces in width. These two properties of 
the energy dispersion become overpowering effects on the 
s-wave gap function's magnitude for larger values of 8, 
despite t' becoming more negative. This is an expected 
density-of-states effect, already pointed out by several 
authors @|§. 

The d-wave magnitudes of Fig. ^(b) are the most strik- 
ing results of all. The magnitudes are highest closer to 
half-filling, again as expected from the Cooper problem 
results, where 8 has values more representative of actual 
physical superconducting systems. For a given value of t', 
the highest gap function magnitude seems to occur at a 
hole density representative of optimal doping for the par- 
ticular superconducting material being modeled by that 
t' value. 

For example, the highest gap function magnitude for 
t' = —0.15 t occurs at 8 = 0.15, which is close to optimal 
doping for I^-zSr^CuO^, whereas for t' = —0.45 t, the 
highest magnitude lies at 8 = 0.4, close to optimal doping 
for YBaaCuaOr-*. 

Further mean field calculations were made with re- 
gard to anisotropy. The previous results were all made 
with tetragonal symmetry, where the relative phase be- 
tween the s- and d-wave components of the solutions 
to the gap equation was found to be tt/2 and when 
the starting matrix contained complex and real ran- 
dom values, respectively. The ground state energies for 
both cases were found to be degenerate, at least for 
the parameter ranges used in these calculations. When 
an anisotropic energy dispersion was assumed, namely 
e(k) = — 2t a cosk x — 2t{,cosk y — At' cos k x cos k y , where 
t a =/= tf,, this relative phase between the s- and d-wave 
components was only found to be 0, even when an ini- 
tial matrix of random complex values was used. As 
the anisotropy increased, ip s increased and ip$ decreased, 
however, both effects were not profound. 

A systematic study of the extended t-J model has been 
presented in the last three sections, beginning with a sim- 
plistic single-pair model in the dilute limit, followed by 
evaluation of a single-pair in the Cooper problem, and 
then ending with a mean field approach. At each step 
along this evolution, we have shown that s- and d-wave 
phases do exist (and indeed co-exist in a mixed fash- 
ion), dependent upon the shape and size of the Fermi 
surface as described by the parameters t' and 8. Varia- 
tional Monte Carlo calculations have already shown the 
presence of a mixed s-d symmetry in the t-J model for 
the case of t' = [^6|. A cross-over was seen where the 
ground state symmetry was pure d-wave for 8 < 0.12 
and an admixture of s- and d-wave for 8 > 0.12 up to 
around 8 = 0.4. Further unpublished calculations show 
that the addition of the t'-term decreases and eliminates 
this mixed phase as t' increases in magnitude. The trends 
derived from the dilute limit, the Cooper problem, and 
the mean field approach as treated in the t-J model have 
now yielded physical reasons for the results seen in the 
VMC calculations. 
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Support for the applicability of the extended t-J model 
in describing the high-temperature superconductors must 
far more importantly come from experimental evidence 
for a gap symmetry dependent upon Fermi surface shape 
and doping. A s urvey of various experimental work is 
made in Sec. VII where the case is made for identifying 
the extended t-J model as a detailed material-specific 
description of superconductivity in the cuprates. 



VII. EXPERIMENTAL SURVEY 

The calculations reported in this paper identify two 
trends for the ground state of high-T c superconductors. 
There is dependence of the form of the gap function on 
the dispersion relation quantified by the second-neighbor 
hopping. There is dependence on the hole concentration. 
These are parameters which vary considerably from sys- 
tem to system. The theory therefore indicates that not 
all high-T c superconductors are alike. Experiments on 
magnetic properties also indicate differences between the 
different systems; the incommensurability of the short- 
range order in the L^-zSr^CuC^ system does not ap- 
pear to be present in the YBa2Cu307 system. Similarly, 
these properties are doping-dependent; the various sig- 
natures of the spin gap are present in underdoped, but 
not optimally doped YBa 2 Cu 3 07^,5. (For a review of 
magnetic properties, see Ref. 35 1.) This section is de- 
voted to a survey of different systems to see if there is 
experimental support for this point of view. Indeed, a 
growth of theoretical work has resulted from these recent 
experimental results, the details of which we present in 
Sec. VII. According to the above results, we should in- 
vestigate each separate high-T c for d-wave, s-wave, and 
s-d mixed wave |3^j3^]. The last possibility requires 
further discussion, s-d mixing can occur in two forms: 
s+id and s+d, depending upon whether the proportion- 
ality constant between them does or does not have an 
imaginary part, respectively. The calculations above do 
not discriminate between these two possibilities. For re- 
cent work on the mixing hypothesis in various models, 
see Refs. p8| , ^9| . Variational Monte Carlo calculations, 
which could in principle show an energy difference in the 
ground state as a function of the phase difference, were 
inconclusive [jL6| . In the s+id state there are no nodes in 
the gap, and time reversal symmetry is broken. Joseph- 
son interference experiments put very strong upper lim- 
its on the admixture of s-wave in any s+id state p0| . 
In addition, the s+id state would lead to domain effects 
which are not observed. (By contrast, the negative re- 
sults in optical dichroism experiments do not test for the 
s+id state 0.) 

In the following discussion, however, we shall concen- 
trate on the s+d state as the most likely mixed pairing 
state, and all references to mixed-wave pairing below are 
to this state only. In this state the d-wave gap nodes, 
located along the k x = ±k y directions, move toward the 



crystal axis k x = or k y — as the relative weight of 
the s-component increases. Which of the two directions 
is chosen depends on the sign of the coupling to the or- 
thorhombic distortion which exists in most high-Tc ma- 
terials. If there is no such distortion, then s + d and s — d 
are degenerate. If a critical size of the s-wave component 
is reached, then the nodes move together and they anni- 
hilate one another. Beyond that point, the system has a 
"hard" gap; there are no quasiparticle excitations. If the 
s-wave component predominates, the d-wave admixture 
may be thought of as contributing some orthorhombic 
distortion to a gap with otherwise tetragonal symmetry. 
Further discussion of the mixed state is found in Ref. [|l6| . 

The purpose of this section is to survey what is known 
about the gap symmetry of different high-T c materials, 
without making the assumption that the gap symmetry 
must be the same in all. Previous surveys have often 
lumped all materials together. This would defeat our 
purpose of looking for material-specific trends. A more 
complete survey, but from a different point of view, may 
again be found in Ref. (ll| . 

The best-studied compound is of course YBa2Cu307. 
It appears from temperature-dependent penetration 
depth measurements A(T) to be quite clear that there 
are low energy quasiparticle excitations in this system 
0. This shows that any admixture of s-wave must be 
less than the critical value. Anisotropy in A suggests that 
pure d-wave is not the gap symmetry in this material, but 
this could also be due simply to the orthorhombic char- 
acter of the crystal structure. 

The Josephson interference experiments of the Illinois 
and ETH groups measure the difference in the phase of 
electrons passing through perpendicular faces of a crystal 
HQ. It is rather accurately measured to be rr, and this 
is in agreement with independent measurements. This is 
likely to reflect very well the phase shift under a 90° ro- 
tation in momentum space. It is important to point out 
that this is consistent with pure d-wave or with the s+d 
state, both of which give n for this quantity. These exper- 
iments rule out the s+id state, which has a shift of 7r/2 
if the two components have equal weight EfJ. The IBM 
tricrystal experiments Q further show that the nodes in 
the gap cannot be more than 30° from the k x = ±k y 
points. This puts constraints on the weight of any ad- 
mixture of s-wave. c-axis tunneling experiments, on the 
other hand, would suggest that the s-wave component is 
not negligible 0]. Indeed, if the critical current along 
this direction is related to the size of the s-wave gap, the 
s-wave part is perhaps 10% to 20% of the d-wave part. 
Nuclear magnetic resonance experiments on this mate- 
rial are inconsistent with a pure s-wave state, and have 
often been interpreted in terms of d-wave superconduc- 
tivity [jl3). The s+d interpretation has not been done; 
however, the fact that anisotropic s-wave [Q and s+id 
p6[ theories can account for the data equally well sug- 
gests strongly that the observations are consistent with 
s+d as well. Temperature-dependent Raman scattering 
on YBa2Cu3C>7 shows most sensitivity in the B\ g ge- 
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ometry [|45| . This supports d-wave pairing Q,[47|] . Some 
limitation can be placed on the proportion of s-wave from 
these experiments, but how much is not clear. Quasipar- 
ticle (Giaever) tunneling experiments performed by point 
contacts [Q , tend to show the presence of states near the 
chemical potential, consistent with d-wave or s+d, but 
it has generally proven to be difficult to unambiguously 
compare theory and experiment (and even experiment 
and experiment) in this area. It appears that quasipar- 
ticle tunneling performed in conjunction with scanning 
tunneling microscopy holds out the best hope of gap de- 
termination, as surface quality and site specificity are 
important issues |49| . At present, however, this method 
gives ambiguous results for the density of states in the 
superconducting state of YBa 2 Cu3C>7 Q. 

Using Ginzburg-Landau theory, Berlinsky and Kallin 
studied the Abrikosov vortex lattice and found that the 
lattice changes from a triangular to an oblique to a square 
spacing with increasing field and sd-mixing parameter 
pl[ . This may explain neutron scattering experiments in 
YBa2Cu307 [p2[ . The interesting point about this work 
from the point of view of this paper is that the s-wave 
component of the order parameter is zero in the absence 
of current flow. Instead, it is induced by the supercurrent 
which is itself a response to the applied field Q . Thus, 
there does not need to be an actual s-wave instability in 
the (translationally invariant) ground state for the effects 
of s-d mixing to be experimentally observable. It may 
be sufficient that there exists a phase boundary between 
pure d wave and and s-d mixed phase and that some 
systems are close to it. 

The second material in terms of detailed gap studies is 
Bi2Sr2CaCu208+a;. Its particular interest is that the gap 
is visible in angle-resolved photoemission. Although this 
experiment measures only the absolute magnitude of the 
gap function, it has an angular resolution not available 
in the Josephson tunneling experiments. Furthermore, 
x can be varied, and this is of pivotal importance from 
our point of view, as the prediction that the s-wave com- 
ponent can develop as the hole doping (a;) increases can 
be tested. The optimally doped compound shows a well- 
developed gap A x along at least one of the crystal axes 
(say, k x = 0) , and a very small or zero gap A xy along the 
diagonal (say, k x = ±k y ) This is in contrast to results 



tion at lower temperatures, and that the anisotropy ratio 
A x /A xy is therefore temperature-dependent. This has 
been successfully analyzed under the hypothesis of s+d 
and is inconsistent with pure d- wave H|. Raman spec- 
troscopy in this compound also supports d-wave pairing, 
again with no clear limitation on s-d mixing ]55| . Giaever 
tunneling experiments are not always consistent with one 
another in Bi2Sr2CaCu20s+2; materials, but again often 
suggest states near the chemical potential. Recent scan- 
ning tunneling microscopy experiments have been seen 
as indicative of an s+d state 56 , a d state |57|, and 



in overdoped Bi2Sr2CaCu208+a; (T c = 75K) 
is found that there is a gap along the k x — 
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Experiments probing gap structure are somewhat 
harder to come by in other materials. Nd2-zCe a; Cu04 
has been the subject of experiments measuring the tem- 
perature dependence of the penetration depth B] and the 
tunneling spectra |59| . Both are consistent with a fully- 
developed gap of s-wave type. Neutron scattering exper- 
iments in La2- a; Sr 2 ;Cu04 detect a gap in the k x = k y 
direction ]60| . This experiment is consistent with s-wave 
or s+d, though the somewhat isotropic suppression of 
scattering would favor pure s-wave, probably of a gapless 
nature. Impurity scattering in a d-wave superconductor 
can also produce rather isotropic suppression ]6l[ |. The 
tricrystal experiment has recently been performed on a 
Tl-based film |32| . It also indicates that nodes in the gap 
of this material must be within 30° of the diagonals in 
the Brillouin zone. 

We summarize this information in Table [|. Also given 
in the table are values of t', the second neighbor hopping 
element, taken from various sources J63|-p5|. The theory 
would generally predict that s-wave should be more likely 
as one moves to the right in the table, as t' decreases in 
magnitude. The doping is not so easy to estimate, but 
it is certain that the overdoped Bi2Sr2CaCu208+:E has 
more holes than the optimally doped Bi2Sr2CaCu20s+a;, 
and therefore should lie in the table as shown. In general, 
the trend is indeed toward more s-wave behavior on the 
right hand side of the table, in qualitative agreement with 
theory. 

A common experimental objection to the s-d mixed 
state is that it involves two superconducting transitions, 
but only one is observed. This subject is treated rather 
thoroughly in Ref. Jl6| and we summarize their conclu- 
sions here. Symmetry analysis shows that there are two 
transitions only if the material is tetragonal. In the or- 
thorhombic case, s-wave and d-wave have the same trans- 
formation properties under the operations of the point 
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group and only one transition will occur. In UPta, the 
only material in which two superconducting transitions 
are known to occur as a function of temperature, they 
are visible in specific heat measurements and in ultra- 
sonic absorption and velocity measurements. In high-T c 
systems, the transition in specific heat has been clearly 
seen only in YBa2Cu307, which is quite orthorhombic. 
In the nearly-tetragonal Bi-based and Tl-based systems, 
the appropriate measurements have not been done. 



VIII. CONCLUSION 

The calculations presented in this paper use the sim- 
ple mean-field method for the investigation of supercon- 
ductivity. There is, in practice, little choice of methods 
if the Hamiltonian has a form which is flexible enough 
to show material-specific behavior. (See, for example, 
Eq. (|7§|).) More sophisticated methods, such as slave- 
boson methods with gauge fluctuations, have been devel- 
oped for the investigation of t- J-type models, but these 
are quite difficult to handle even for the simplest mod- 
els. Finite-size calculations are limited to lattices whose 
linear dimensions are comparable to, or smaller than, a 
coherence length. Quantum Monte Carlo (QMC) calcula- 
tions (again on simpler models than the current one) have 
severe sign problems in the regimes of interest. VMC 
calculations may be possible, but this is left for future 
investigation. To the extent that we have been able to 
compare our results with the VMC results, there is qual- 
itative agreement. Finally, our calculations in the spin- 
interaction model are cruder than the strong-coupling 
calculations which have been done in the spin-fluctuation 
model. The spin-interaction model is genuinely micro- 
scopic, however, not semiphenomenological. Compensat- 
ing for the disadvantage of the approximate treatment of 
interactions is the advantage that the method is physi- 
cally transparent - the reasons for the trends which ap- 
pear are all evident, and much of the paper has been 
devoted to the explanation of these trends. This trans- 
parency is not always a feature of highly numerical meth- 
ods such as QMC or VMC. 

Thus the trends towards d-wave with increasing t' can 
be seen as softening angular variations in the pair wave- 
function, whereas the trend towards s-wave with increas- 
ing doping is seen as due to a relaxation of phase-space 
restrictions. The trend towards higher optimal doping 
with increasing t' is a density-of- stated effect. There is 
evidence for all of these effects in experiments. The most 
dramatic of these is the change from strong d-wave in 
YBa2Cu307 to s-wave in Nd2- a; Ce 2 ;Cu04. 

This analysis allows us to evaluate qualitatively the ac- 
tual relevance of the calculations to real materials. For 
example, let the ideal Fermi surface of the Cooper prob- 
lem be replaced by the actual momentum distribution. 
The latter will be much flatter, with a relatively small 
discontinuity at the Fermi surface (assuming a Fermi 



liquid normal state) or no discontinuity (Luttinger liq- 
uid normal state). Therefore the phase space restric- 
tions will be weaker in the real system. This would fa- 
vor s-wave and move the crossing point of s-wave and 
d-wave at n = 0.55 (at t' = 0) in Fig. | to higher val- 
ues with increasing t' . Relaxation of the very large U 
(strong short range Coulomb interaction) approximation 
will also favor s-wave, as the constraint of absolutely no 
double-occupancy has no effect on <i-wave, whereas it 
suppresses s-wave pairing. By the same token, the ne- 
glect of the long range Coulomb interaction in our model 
probably prejudices our results toward d-w&ve. The usual 
frequency-mismatch arguments so important for super- 
conductivity with electron-phonon interaction do not ap- 
ply here, since all interactions are instantaneous. It is 
precisely for that reason that our results fly in the face of 
conventional wisdom that magnetic mechanisms can lead 
only to d-wave pairing. In the absence of a cutoff in the 
interaction, all of momentum space comes into play, and 
the kinetic energy of the pair plays an important role. 
This brings s-wave superconductivity into the picture. 

The trends which arise from our simple analysis of this 
flexible model appear to describe some of the otherwise 
puzzling and contradictory results on different high-T c 
materials. They show that the ground state of different 
materials may be different in the context of a magnetic 
mechanism. 
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